Electronic transport driven spin-dynamics 
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We propose a model to explore the dynamics of spin-systems coupled by exchange interaction to 
the conduction band electrons of a semiconductor material that forms the channel in a ferromag- 
net/semiconductor/ferromagnet spin- valve structure. We show that recent observation of the novel 
transient transport signature in a MnAs/GaAs/MnAs spin- valve structure with paramagnetic Mn 
impurities [D. Saha et a/., Phys. Rev. Lett., 100, 196603 (2008)] can be quantitatively understood 
in terms of current driven dynamical polarization of Mn spins. Using our model of spin polarized 
transport through Schottky barriers at the two ferromagnet /semiconductor junctions in a spin- valve 
structure and a dynamical equation describing the paramagnetic impurities coupled to conduction 
band electrons we explain the scaling behaviour of observed transient features such as the magnitude 
and time-scale with temperature. 

PACS numbers: 



I. INTRODUCTION 

Spin-dynamics are finding wide interest in recent years 
as they usher in the prospect of finding their applications 
in quantum computing p], 0], data storage [3[, biolog- 
ical and medical applications 0, [Bj], etc. The success 
of these goals are contingent on being able to control 
and detect the spin degrees of freedom concerned. The 
succesful experimental demonstration of spin-injection 
from ferrmoagnetic contacts into semiconductors 0, 0, 
and of controlling the spin-injection and extraction by 
spin-valve structures in semiconductors Q enables the 
possible control of surrounding spin-degrees of freedom 
by spin-polarized carriers inside the semiconductor stru- 
tures. Understanding the ensuing spin dynamics by 
the spin-polarized carrier is of immense importance to 
all novel electrically driven spin-based device concepts 
ranging from devices employing non-interacting spins Q 
to those employing strongly interacting spins in mag- 
nets [To(] . 

We have developed a model to treat the non- 
equilibrium non-interacting spin dynamics driven by 
spin-polarized electrons injected into and extracted from 
a semiconductor channel via Schottky barriers between 
the semiconductor and ferromagnetic contacts arranged 
in a spin-valve configuration. By modeling the transport 
through these two Schottky barriers and solving them 
together in conjunction with a minimal transport model 
for the semiconductor, we can estimate the spin- flip cur- 
rent inside the semiconductor channel and thereby calcu- 
late the ensuing dynamical evolution of non-interacting 
spin-systems having exchange-interaction with the chan- 
nel electrons. 

We apply our model to analyze recent experimental ob- 
servations [ll| in MnAs/GaAs/MnAs lateral spin- valve 
structure with paramagnetic Mn impurities embedded in 
GaAs (fig. [J). We attribute the origin of novel transient 
features such as the decay in the anti-parallel current 
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FIG. 1: Schematic structure of the Spin-charge Transducer 
(SCT). Lp( a ), L, W, and, x denotes the length of the polarizer 
(analyzer), separation bewtween contacts, width of the entire 
structure, and, thickness of GaAs channel respectively. The 
magnetization of the analyzer can be reversed back and forth 
to to make it parallel (P) and anti-parallel (AP) to that of 
the polarizer. 



over time (fig. [2]) to dynamical polarization of paramag- 
netic Mn impurities (inset of fig. Ej) and can successfully 
explain the scaling of the transient characteristics such 
as the magnitude and time-constant with temperature 
(fig. E]). In this regard we point to the suppressive role 
played by the other electronic spin-flip mechanisms in 
GaAs, which do not involve Mn impurities (fig. [6|), to 
reduce the effect at elevated temperature. 

Our paper thus serves a two- fold purpose: (i) it identi- 
fies novel transport signatures pertaining to the underly- 
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FIG. 2: Transient variation of terminal current / for an 
applied bias V = 25mV during parallel (Ip'.top) and anti- 
parallel (I ap' bottom) configuration at different temperature 
for an SCT structure with dimensions: L p = 5/xm, L = 1/mi, 
L a — 30/xm, W — 50/xm, and, x — 0.5/xm. 



ing physical process of dynamical polarization of param- 
agnetic impurities embedded in a material carrying spin- 
polarized electrons which are coupled by exchange inter- 
action to the impurities, and (ii) it identifies the effect of 
temperature and other electronic spin-relaxation mecha- 
nisms on the mentioned signatures. 



II. THEORY AND CALCULATION 

A. Spin transport through the SCT 

The basic electronic transport mechanism through the 
SCT is that electrons tunnel through two back-to-back 
Schottky barriers at the two MnAs/GaAs junctions, 
which are connected via a diffusive semiconductor chan- 
nel (fig. [3j). To model the tunneling Schottky barrier 
the magnetized MnAs contacts the model equations in 
refs. pjl, needs to be modifeid to consider that (i) the 
magnetized MnAs density of states and Fermi velocity for 
the up and the down electrons are different around the 
chemical potential, and (ii) the up-spin and the down- 
spin electrons might see different electrochemical poten- 
tials Ktp(a) an d Vdp(a) respectively while tunneling across 
the forward (reversed) biased polarizer (analyzer) junc- 
tion but the same electrostatic potential V esp ^ as the 
electrostatic potential depends on the total charge den- 
sity. We introduce a contact polarization parameter P c 
to consider the first effect which implies that even if both 
the up-spin and down-spin electrons see the same elec- 
trochemical potential while tunneling across the barrier, 
the ratio of majority to minority spin-component of the 
current is, (1 + P c ) : (1 — P c ). In all our calculations we 
will follow the convention that when the two magnetic 
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FIG. 3: (a) Equillibrium band diagram: the GaAs channel is 
coupled to the MnAs contacts through Schottky barriers. Ef 
and E c are the equillibrium chemical potential and conduction 
band edge of the GaAs respectively, (b) Band diagram under 
bias: The polarizer Schottky barrier is forward biased and the 
Analyzer Schottky barrier is reverse biased. The electrons are 
being injected into (extracted from) the GaAs channel, main- 
taining a chemical potential \i c han-> from (to) the polarizer 
(analyzer) MnAs conatact having chemical potential fi p (/x a ) 
(c) Equivalent circuit diagram for SCT with non-magnetized 
contacts: the applied voltage V drives a terminal current / 
through two diodes connected in series by a resistor. 

MnAs contacts have parallel magnetization the majority 
(will be denoted by M) and minority (will be denoted 
by m) spin-components of both the contacts are up-spin 
(will be denoted by u) and down-spin respectively (de- 
noted by d) . Then it follows that when the contacts have 
anti-parallel magnetization the majority and the minor- 
ity spin-components of the analyzer will be the oppo- 
site of that of the parallel case (i.e. down-spin and up- 
spin respectively) as it's magnetization is reversed during 
the experiment. The second effect will directly modify 
the tunnel current expressions in refs. pjj 05 in such 
a way that the electrostatic potential essentially enters 
the terms that comes from the shape of the barrier and 
the two electrochemical potentials enters the terms that 
involves the position of the Fermi levels. After making 
these modifications we finally get the expressions for ma- 
jority and minority spin-component of current through 
the polarizer (analyzer) I p ( a )M an d I p (a)m respectively: 

RB ( F B} 

Ip(a)M = 

0.5(1 +P c )L p{a) WJ2h K \Ves P {a),V Mp{a) ){l) 
Ip(a)m = 0-5(1 - P c )L p ^WJ^ c f ^ FB \V esp ( a ),V rnp ( a ))(2) 

where W is the width of the n-GaAs channel (fig. [1]) 
and £p( a ) is the length of the polarizer (analyzer) and 
(M, m) = (^, d) for both the contacts during parallel 
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magnetization and (M, m) = (d, only for the analyzer 
during anti-parallel magnetization. The current density 
is given by: 

_ b FB(RB) 
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where, ?/ G {M, m}, is the Boltzman constant, g is 
the electronic charge, and T is the temperature. Now for 
forward bias: 
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and for reverse bias: 
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Here, £ o = 2q^N D1 /2e/a with a = 2V2m*/^- ft is 
the Plank's constant, m* is the effective mass of GaAs 
conduction band electrons, e is the dielectric permittiv- 
ity of GaAs, Njji is the doping density of the n+GaAs 
region (figs. [Hand 03(a)), Ef c i is the position of the equ- 
ilibrium Fermi energy of the n+GaAs above it's conduc- 
tion band-edge, and (j) sm is the Schottky barrier height at 
equillibrium, which will be used as a fitting parameter. 
In the expressions above all the electrochemical and elec- 
trostatic potential are not completely independent due 
to the restriction of charge neutrality. Assuming that 
the difference between the up-spin and down-spin charge 
density to be smaller enough than the total charge den- 
sity we can write: 



V e3 a(p) = (VMa(p) + Krm(p))/ 2 



(8) 



i(p) ~r v ma{p) 

Finally, electrostatic potentials across the polarizer, an- 
alyzer and the channel (V c han) should add up to applied 
voltage V and the currents through the polarizer (J p ), 
analyzer (I a ) and channel (I c han) should all be equal to 
the current in the external circuit: 
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FIG. 4: Equivalent circuit diagram for the SCT with magne- 
tized contacts 



where, x, Nr)2, and v are the thickness, doping density 
and mobility of the nGaAs channel respectively (figs. [1]) 
and L is the separation between the contacts. 

From equations[4l[6]and[7]we see that the electrochem- 
ical potential across the reverse biased junction only en- 

_ FB(RB) , v | 

ters through {1 — e Cl 91 y y<p) 1 } which becomes 1 for 
the reverse biased polarizer voltage of interest. So we can 
assume that for the reverse biased junction (i.e. polar- 
izer) the exact value of the electrochemical potential V yp 
is of no consequence and will always be set equal to elec- 



trostatic potential V ei 
i.e. 

Vmp 



in all the calculations to follow, 
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B. Estimating the spin-flip current 

As the majority and minority electron spins enter and 
flow through the channel material, different spin-flip 
mechanisms can transfer some electrons from the up-spin 
stream to the down-spin stream and vice-verca. As a re- 
sult, the extracted up(down)-spin current I au u) might 
not be the same as the injected up(down)-spin current 
Ipu(d)- The difference is what we will call the spin-flip 
current: 



Lpd 



(15) 



A schematic circuit diagram with different current com- 
ponent discussed so far is shown in Fig. 0J 

Parallel configuration: We will start by analyzing the 
data for the parallel configuration. In this configuration 



we can rewrite the eqs. [T] and |5J 

I pu = 0.5(1 + P c )L p Wjn%{V esp , V up = V esp ) (16) 

I pd = 0.5(1 - P c )L p WJ*%(V esp , V dp = V esp ) (IT) 

I au = 0.5(1 + P c )L a WJf£(V esa ,V ua ) (18) 

I ad = 0.5(1 - P c )L a WJ™(V esa ,V da ) (19) 

In all following we will not consider the negligible time- 
dependence of parallel current data and will assume that 
the it has a constant value over time. Then it follows 
that there is no significant spin-flip dynamics taking place 
during parallel configuration and we can assume spin-flip 
current to be zero. For zero spin-flip current we get from 
eq. [13 [Hand CCD 

¥ = ¥ = ^ 

-Lpd -'-ad -L -*c 

Now from eqs. [8] [H [J9] and [20] we get: 

V ua = Vda = V esa (21) 

Substituting [2T1 and fTGl - [T9l into fTUl we find that the par- 
allel current: 

I P = L p WJ*%(V eap , V esp ) = L a WJ? c %(V eaa , V esa ) (22) 

is independent ofP c and only parametrically depends on 
</> sm . For a particular value of parallel current 7 and 
a particular value of (j) sm we can solve eq. [22] to deter- 
mine V esp and V esa . From the knowledge of I c han = I 
we can determine V c han from eq. [13] Subsequently, we 
can calculate V from eq. [9] for certain values of 7 and 
(j) sm . Conversely, by aplying iterative solution technique 
on these steps we can determine (j) srn from the experimen- 
tal value of 7 and V for parallel case and the calculated 
value is (j) srn ^0.6 which is very close to the known value 
in literature [15[ and fits the data for all the experimental 
temperatures. 

Antiparallel configuration: The prominent transient be- 
haviour of the AP current indicates significant spin-flip 
dynamics. With the value of (j) sm estimated from the P 
configuration data we will now proceed to estimate the 
spin-flip current variation over time. We first re- write 
eqs. [T]and [2] for the AP case: 

I pu = 0.5(1 + P c )L p WJfi(V esp , V up = V esp ) (23) 
I pd = 0.5(1 - P c )L p WJ*i(V esp , V dp = V esp ) (24) 

hu = 0.5(1 - P c )L a Wjf£ {V eaa , Ka) (25) 

hd = 0.5(1 + P c )L a WJ™(y esa ,V da ) (26) 
By adding eqs. [23l and [2^1 we get the AP terminal current: 

lAP = Ipu + Ipd ~ LpW J^ch(Vesp, V esp ) (27) 

Now for a particular value of AP terminal current Iap we 
determine the value of electrostatic potential V esp using 
eq. [27] Then for a particular value of contact polarization 
P c we determine the up-spin and down-spin components 
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FIG. 5: Time variation of spin- flip current, I s f at different 
temperaures during antiparallel configuration, inset: Varia- 
tion of spin-flip current I s f with antiparallel terminal current 
normalized by parallel terminal current (Iap/Ip) at those 
temperatures. All the results are for same the voltage and 
dimension values as that of fig. [2] and P c = 0.49 18]. The 
values of the other parameters used for this calculation are: 
N D1 = 1 x 10 19 cm- 3 , E fcl = 400meV, sm « 0.6eV [3, 
and, N D2 = 8.6 x 10 16 cm- 3 . 

of the current (I pu and I p d respectively) injected through 
the polarizer terminal from: 

7^ = 0.5(1 + ^)7 (28) 
Ipd = 0.5(1 - P C )I (29) 

We next determine the two spin-components 7 an and I a d 
of analyzer current from[15l At this stage we know the 
values of (j) sm from parallel data, the currents I au and 
I a d for a particular value of P c and solve three eqs. [H [251 
and f26l for three unknown variables V esa , V ua and V pa . 
Finally we determine V c han from eq. [13] and [10] With 
different potentials V esp , V esa and V c han known we readily 
determine the applied voltage V from eq. [9] By execution 
of all the steps mentioned above we finally deduce the 
value of V for arbitrary values of P c , 7 s j, and 7 using the 
value of (j) srn determined from the parallel data. Again, 
by iterative solution technique we can use these steps to 
determine the value of I s f for a particular value of 7, V, 
and P c . The results of such a calculation are shown in the 
inset of fig. [5l Using these results (in the inset of fig. [5]) 
and the values of the experimental antiparallel current 
(fig. [2]) we determine the time variation of I s f which is 
reported in fig. [5] 

C. Impurity dynamics 

The Mn impurities in GaAs are compensated. Hence, 
they have a single negative charge and are in S = 5/2 
spin state [l6[ . We will estimate the impurity spin polar- 
ization in terms of the average z-component value (S z ) 
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of impurity spins: 

(S z ) = 2^+5/2 + 2^+3/2 + 2 Fl/2 

where F n is the probability of an impurity atom being 
in S z = n state. When the impurity-spins are exchange 
coupled to the electron-spins the ensuing electron spin- 
flip current lfj n can be written as (by extending the 
equations in [17[ to S=5/2 system): 

If s n = l^J 2 Nj J D u (E)D d (E) [(-F-5/2 +-F-3/2+ 
F_i/2 + F +1/2 + F +3/2 ) f u (E){\ - f d (E)}- 

{F-3/2 + F-l/2 + F+1/2 + -F+3/2 + F +5 / 2 ) x 

f d (E){l-f u (E)}]dE (31) 



where, D u ^ d ) and f u ^) are the density of states (DOS) 
and average energy distribution functions of up(down)- 
spin electrons inside the channel, Ni is the total no. of 
impurities, and, h = 2ttH. As the impurities get spin- 
polarized their dynamics can be described as: 
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Here, 7/ is the relaxation rate of the impurity-spins in- 
dependent of electrons and 
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where, 

27T /* 

A =Y 2 J j2jD «^/«(i - / d )d£ 

27T /* 

S = —2 y J 2 D u D d f d {\ - f u )dE 
Now from eqs. [30j [31j and, [32] we can write: 

J_,. r(t) ^ +7l ( Sl(()) (34 ) 

As the spin-lattice relaxation of Mn impurity spins are 
negligible within the experimental temparature range [l6[ 
7/ « 0. Consequently, it follows that I^j n (t — > 00) = 
i.e. the spin-flip current component involing Mn impu- 
rities have zero contribution to the steady state spin-flip 
current. It is the other component of the spin- flip current 
which involves conduction band electron spin- flip due 
to the spin-orbit interaction of GaAs, that gives rise to 
steady state spin-flip current: 1*^ = I s f(t — > 00) and is 
always flowing as a 'background' spin- flip current. There- 
fore we can find the spin- flip current due to Mn impurities 
from the results of fig. [5] as lf f n (t) = I af (t) - I sf (t 00) 
and the results are shown in fig. [6] 

The GaAs channel being diffusive, we can assume that 
the energy distribution function of the up (down) elec- 
trons f u u) t° be Fermi functions with temperature T 



and chemical potentials fi u (d)- We assume that \i u — fid 
is small and does not vary significantly over time. Then 
we have, due to charge neutrality condition inside the 
channel, the constraint (fi u + /id)/2 = Ef C 2, where Ef C 2 
is the Fermi level of the n-GaAs channel correspond- 
ing to the doping of the channel (and not to be con- 
fused with Ef c i which is the Fermi level of n+GaAs re- 
gion). Also the GaAs having spin-degenerate bands (i.e. 
D U (E) = Dd(E)) we can consider the density of states 
inside the channel to be uniform over the small energy 
window and define D = D u (Ef C 2) = Dd(Ef C 2). With 
these simplifications we can calculate the time variation 
of I^j n and (S z ) from eqs.[30l[3TJ and, [32] for a particular 
value of JD and fi u — ^d with the values of Ef C 2 and T set 
by the experimental conditions. Such a calculation was 
performed by varying JD and fi u — fid to match I^j n cal- 
culated from fig. [5] for T = 10K and the results are shown 
in fig. [6] As the temperature is raised the value of JD 
should, in principle, remain unchanged and so a similar 
calculation was performed by varying only fi u — fid and 
keeping the value of JD unchanged to the previously de- 
termined value to successfully match the I^j n calculated 
from fig. [5] for T = 15K. These results are also shown in 

fig- El 
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spin and down-spin electrons (fi u — fid in fig. [6]) by deter- 
mining the spin- flip conductance g so = I^/(fi u — fid) and 
then equating the conductance value to the expression 
^-DWx(L + L p + L a )^r-. The estimated spin-relaxation 
time Ts of the electrons independent of Mn impurities is 
found to be 2ns (Ins) for T=10K (15K), which are very 
close to to the reported values [~L9^ . 

Furthermore, the strength of the exchange interaction 
(J/ag) between conduction band electrons and the Mn 
impurities is estimated to be in the order of leV (from 
the value of JD in fig. [6] with the value of D calculated 
using Ef C 2 = lOmeV corresponding to channel doping) 
which is an order of magnitude larger than the estimated 
s-d exchange interaction at band-edge in I20I1 and two 
orders of magnitude larger than the same in [2l| . 



FIG. 6: Transient variation of I s f n calculated from the impu- 
rity dynamics (solid and dashed lines) and from the calcula- 
tion corresponding to fig. [5] (circles and triangles) (the results 
corresponsing to 10K has been shifted upwards by 2/nA. In- 
set: Transient variation of the average z-cmponent (S z ) of Mn 
impurity spin at different temperatures. The parameters val- 
ues for calculation corresponding to the impurity dynamics 
are: JD = 4.8 x 10" 3 and fi u - fi d = 0.25mV (0.125mV) for 
T= 10K (15K). 



III. DISCUSSION 

Comparing the results in fig. [6] at different temper- 
atures we find that the difference between the chemical 
potentials of the up-spin and down-spin electrons fi u — fid 
decreases with increasing temperature. This is to be ex- 
pected as the spin-orbit component of the spin-flip cur- 
rent I* J due to electron spin-relaxation mechanism in 
GaAs is known to increase with increasing temperature 
from magnetoresistance measurements on spin- valves [?} . 
Such increase in Jjj 'loads' the chemical potential differ- 
ence between up-spin and down-spin electrons fi u — fid 
that acts as a source to this component. As fi u — fid de- 
creases with increasing temperature the initial value of 
I^f n , which is proportional to / (f u - fd)dE ~ (fi u - fi d ) 
(eq. IA2p . is expected to decrease according to our theo- 
retical analysis and is also observed in the experiment. 
On the other hand, it can be shown that the time-scale 
of the transient variation, according eq. [32J, is inversely 
proportional to $ {f u {^ ~ fd) + fd{^ ~ fu)}dE which can 
be shown to be equal to 2ksT for such small value of 
fiu — fid- Hence according to our theory we expect the 
time-scale to be inversely proportional to the tempera- 
ture, which is observed in the experiment: the value of 
time constant changed from 120ns to 80ns (their ratio 
being 3 : 2) as the temperature was changed from 10K 
to 15K (their ratio being 2:3). 

We estimate the spin-relaxation time of the GaAs con- 
duction band electrons, independent of Mn impurities, 
from the steady state spin-flip current values in fig. [5] 
(J|j) and the chemical potential difference between up- 



IV. CONCLUSION 

We have studied electronic transport driven non- 
equilibrium spin-dynamics arising from exchange interac- 
tion between the conduction band electrons and param- 
agnetic Mn spin impurities embedded in the GaAs chan- 
nel of a MnAs/GaAs/MnAs spin- valve. We believe that 
we have sufficiently and uniquely elucidated the underly- 
ing physical processes leading to the experimental obser- 
vations made in Ref. as we can match several experi- 
mental features (parallel current values at 3 temperature 
settings, initial values and time-constants of I^j n (t) at 
2 temperature settings extracted from the experiment) 
using a few fitting parameters (0 sm , JD and 2 values of 
— /id at 2 different temperatures). From our calcu- 
lations we argue that spin polarized carriers can signif- 
icantly polarize the surrounding spin-systems and that 
the dynamical polarization has an impact on transport 
characteristics. We also observe that the other electronic 
spin-relaxation mechanisms suppresses the mentioned ef- 
fect. In this regard, we would argue that using semicon- 
ductor materials having larger electronic spin-relaxation 
time such as Si, its spin-relaxation time being at least two 
orders of magnitude larger than GaAs [22[, can improve 
the polarization of the impurity spins by reducing the 
'loading' of fi u — ^d given that the spin-injection efficiency 
does not suffer. Nevertheless, the degree of polarization 
acheived by such a small current is quite remarkable if 
we consider the fact that, to attain such polarization one 
would have to otherwise apply a large magnetic field (e.g. 
~ 2 T to get a polarization of ~ 0.8 for an S = 5/2 system 
at 10K). This fact leads us to suggest that spin-polarized 
electronic current can emerge as a good candidate for 
all the experimental investigations and applications that 
necessitates manipulation of various spin-systems. An in- 
teresting extension of this work would be to study the ef- 
fect of transport on weakly-interacting spin-systems such 
as a magnet near it's Curie temperature, which might 
have a prospect for a realization of electrically controlled 
magnetic ordering. 
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APPENDIX A: DERIVATION OF EQ. \M 



to which we apply eqs. l32l and l33l and get: 



Applying the sum rule: 

+5/2 
n=-5/2 

to eq. 02] we get the result: 



I% n = ^2nJ 2 Nj 



, J D u D d {(f u 
F +5/2 f u (l-f d )}dE 



Now from the definition in eq. [ 

T 



d(S z ) 
dt 



+5/2 
+3/2 
+1/2 
-1/2 
-3/2 
-5/2 



d_ 
~dt 



(Al) 



- fd) + F_ 5/2 f d (l - f u ) 
(A2) 

| we can write: 

^+5/2 
^+3/2 
F +l/2 
F-l/2 
F-3/2 



(A - 5) - AF +5/2 + £F_ 5/2 



^2nJ 2 J D u D d {(f u - f d ) + F_ 5/2 f d (l - f u ) 
-F +W2 f u (l-f d )}dE->y I (S z ) 



-5/2 



Substituting eq. IA2I in the above equation we get the 
desired result of eq. EH 
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